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Abstract 

Inelastic scattering of ultra cold neutrons on bubbles with the size of nanometers is considered. 
It is shown that neutron-bubble cross section is large and sensitive to different vibration modes 
of bubbles. This process could be used for study of dynamics of nano-size bubbles, and for new 
methods of ultra cold neutron production if appropriately sized bubbles with sufficient density can 
be created. 
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INTRODUCTION 



The understanding of interactions of ultra cold neutrons (UCN) with materials is very 
important for many applications (see, for example [l| and references therein) including 
fundamental neutron physics, neutron scattering, and condense matter physics . A special 
interest in interactions of UCN with nano size particles as an explanation of the possible 
loss of UCNs during a storage in a trap and as a possible new method for UCN production 
has been emphasized in papers 2|, |3j. A calculation of neutron interactions with large 
particles is a very complicated problem since it requires detailed description of internal 
excitation modes for these particles, which is impossible in general, and can be done only 
numerically for some cases. To avoid these difficulties in general understanding of properties 
of interaction of UCN with nano-particles, we consider a bubble of nano-size in liquid for 
which the neutron scattering problem may be solved analytically. In particular, we consider 
bubbles in liquid helium which is a "natural" environment for UCN applications. 

Let us recall that a characteristic scale of kinetic energy of UCN could be as low as 
at the level of nano-electronvolts. These energies corresponds to neutron wavelengths at 
the nano-meter scale. For example, neutron with the kinetic energy E n = 100 neV has a 
wavelength A = 90.4 nm, and the energy E n = 1 neV corresponds to the neutron wavelength 
A = 904 nm. Then a typical neutron scattering parameter on the nano-size particles is 



(kR) ~ 0.07 — 10, where k = y/2mE n /fi is neutron wave vector and R is size of a target. 
This parameter gives a characteristic range of partial waves required for the cross sections 
calculations. Therefore, in general case UCN cross sections with nano-particles can not be 
calculated in the low energy approximations (kR <C 1), but contributions from many partial 
waves should be taken into account. In other words, very slow neutrons can demonstrate 
"high energy" scattering properties when they interact with nano-particles. 

In this letter, we present results of calculations of inelastic neutron cross sections on the 
bubble of the nano-meter sizes. Some useful results for elastic cross sections are discussed 
briefly since all details for elastic scattering could be found in quantum mechanics textbooks. 
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ELASTIC SCATTERING ON SPHERICAL TARGETS 



For calculation of elastic neutron scattering on the large objects we can use similarity 
between cold neutron scattering and light scattering on transparent objects in optics. In our 
case, it is convenient to use neutron Fermi potential which is related to refractive index in 
optics. The values of Fermi potential for a medium with an atomic density N and neutron 



coherent scattering length 6 is 



v = 2 -^m„, (i) 

m 

where m is a neutron mass. This gives Fermi potential for neutron propagation in liquid A He, 
3 He and solid deuterium as 16 neV, 28 neV and 104 neV jsj, correspondingly. Therefore, 
ultra cold neutrons in helium can feel a potential steps or a potential well of about ~ 20 neV 
(~ 80neV), when they are scattering on bubbles (deuterium spheres), correspondingly. One 
can see that for neutron elastic scattering in this approach there is no differences between 
neutron scattering on bubbles and on spherical solid objects, which is not the case for 
inelastic scattering. To describe this scattering process, one can use a standard results for 
scattering on square- well potential Vq. Then two separate cases with essentially different 
scattering regimes can be considered: a "high energy" regime where parameter kR ~ 1 and 
a "low energy" regime where kR -C 1. The first case results in a total cross section as a sum 
over all important partial waves with phase shifts 5i for angular momenta / < L max : 

L 

k 2 



Att 

a =^ ]T(2/+l)sin 2 ^, (2) 

1=0 



where 



and 
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Here ji(z) = yV / (2z)Ji + i/ 2 (z) and ni(z) = Jtt/ (2z)Ni + \/2(z) are spherical Bessel functions, 



and k = ^2m(E n + Vq)/U. 

In spite of the fact that the potential well is extremely shallow, the size of the well 
can be rather large. Therefore, neutron scattering can show resonance behavior due to a 
sensitivity to a possible virtual (or real) bound states (for neutron bound states in matter 
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FIG. 1: Elastic total cross section of neutrons with energy E n = 100 neV as a function of a bubble 
radius in liquid helium. 
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FIG. 2: Elastic total cross section of neutrons with energy E n = 1 neV as a function of a bubble 
radius in liquid helium. 

sec ) in large bubbles when E n < Vq and 2mVoR 2 /h 2 ~ 1. Numerical calculations 

confirm such phenomena: the cross section as a function of the size of the target is a smoothly 
growing function when neutron energy is much larger than the potential (for example, for 
Vo = 16 neV and E n = 100 neV, see fig. (Q])), but it shows irregular (resonance) behavior 
when neutron energy is much smaller than the potential (for example, E n = 1 neV, see fig. 
(J2])). (For the considered parameters Vo, E n and A ~ R it was required to take into account 
contributions from high orbital momenta up to L max ~ 20.) 

For a "low energy" regime with (kR) 1, one can consider only zero orbital momentum 
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contribution using low energy expansion for the scattering phase. Then, the cross section is 

a ~ ^R 2 (K R) 4 ~ R 6 ~ JV t 2 , (5) 
9 

where i^"o = v^^Vo/ft, and iVt is a total number of nuclei in the target (or the total number 
of nuclei which could fill up the volume of the bubble). One can see a coherent enchantment 
factor N t for elastic scattering cross section. However, the condition required to obtain 
Eq.flHD, (kR) <C 1 is equivalent to the condition A ^> R. Therefore, it does not lead to a 
large value of coherent cross section in compare to a geometrical one (it is suppressed by 
the factor (KqR) 4 as seen from Eq.flS])). 

Since that elastic cross section for scattering of UCN on nano-particles shows regular and 
expected properties, one can conclude that it is not large enough to be considered for the 
mechanisms of UCN energy variations suggested in j^, 3]) where "a sufficiently large cross 



section for coherent interaction" |2j was required 



INELASTIC SCATTERING ON BUBBLES 



Inelastic scattering of UCN on a bubble in a liquid can be calculated using a general 
formula for inelastic cross sections [9J with a transition of the target from an initial i to a 
final / state: 

da fi =^\{fp'\U\ip)\ 2 dtt. (6) 

Assuming that energy of neutrons is high, compared with the neutron binding energy in 
the bubble, one can approximate incoming and outgoing neutron wave functions by plane 
waves. It should be noted that this approximation is always [9] correct if \U\ <C h 2 /mR 2 , 
where R is a radius the bubble. Therefore, it can be safely used for bubbles with radii less 
than tens of nanometers irrespectively to the binding energy condition. Then, the cross 
section can be written as 

2 dn', (7) 

— * — # 

where ipf(r) and ipi(r) are final and initial wave functions of the bubble, and q = k' — k 
is transferred wave vector (k' = p'/h and k = p/%). It is convenient to re- write the cross 
section in terms of q using the identity qdq = (kk' /2ir)dQ' . Then, taking into account that 



da 



fi 



m p 



v 



Ue- iqr i/j* f ^idrdf 



5 



the potential inside the bubble has a constant value Vq, one obtains 

m 2 V r 2 



da 



fi 



'0 



qdq, (8) 



where integration over r is performed inside the volume of the (deformed) bubble. (Since we 
consider neutron propagation in medium, it is natural to assume the reference zero Fermi 
potential for the medium. Therefore the only non-zero effective potential is the potential 
inside the bubble volume.) The shape of the bubble can be described in spherical coordinates 
by equation 1(3] 

R = R + '£a lm Y£ n {e,<f>), (9) 

lm 

where Rq, is a radius of a spherical (non-deformed) bubble, ai m are deformation parameters 
in an expansion by spherical functions Yi m (9, 4>). Assuming small deformations (a>i m <C Rq), 
one can describe dynamics of the bubble in terms of harmonic oscillators in the space of 
deformation parameters [h]]. Then, the integrant for the integration over the space of 
deformation r in the Eq. (jSJ) is proportional to the product the first order of the corresponding 
deformation parameter, and harmonic oscillator functions ipi =n and ipf= n ±i, which results 



in yhn/2Mu for the integral over r. Here n is an integer number, M is a effective mass of 
the bubble, and u is a frequency of the bubble oscillation. 

In this paper, we restrict ourselves by first simple modes of the possible deformations: the 
breathing mode (I = 0) and the surface mode (1 = 2). (The dipole mode with I = 1 describes 
a shift of the bubble in the space and does not lead to a deformation.) For these two modes, 
one can calculate integrals over the r space (/ e~ iqr dr) inside of the deformed bubble with 
the bubble shape described by Eq.([9]). This integral is coupled with harmonic oscillator 
integral by deformation parameters. The harmonic oscillator integral is not equal to zero if 
it is contains linear terms of deformation parameters. Therefore, one can expand the e~ tqr 
function and keep only the part which is proportional to the first order of deformation (linear 
term) in the expansion. Then the matrix element in Eq.® can be written as 



M = I I e-^^drdr = R 2 J J^-^ mode ^), (10) 



2Mu 

where £ = qR®. The function $&(£) for the breathing mode is 



$6(0 = 47^, (11) 
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and the function $ s (£) for the surface mode is 

= ^[S^cose - (3 - e) sin^. (12) 

The Eg. pip is a direct result of integration over bubble volume modified by the first term 
in the radius expansion of Eq.([9]). To obtain Eq. (fT2l) . a simple symmetric ellipsoidal shape 
(022 = 0.2-2 and a 2 \ = a,2~i = 0) of the bubble deformation was assumed. However, the 
assumption about symmetry of deformation is not important for our considerations because 
different shapes may lead only to slightly different coefficients in functions 

Then, differential cross sections for inelastic scattering of UCN on bubbles can be written 

as 

'« = *W (^) "-(Of* « 13 » 

for the case when energy transfers from the neutron to the bubble (down-scattering cross 
section), and as 

m 2 V2Rl (%(n + 1)\ ^ 2 ,,,,,, 

for the case when energy transfers from the bubble to the neutron (up-scattering cross 
section). It should be noted that due to rotational invariance these cross sections are in- 
dependent on m-number in Eq. ([H]). For example, for the case of Eq. fll2l) this independence 
comes automatically after integration over angular parameter 0. The total cross section 
could be obtained by integration of da(C,) from £ m j n until ^ ma x, taking into account that 
k 2 — k' 2 = ±2mu}/h, where "+" at the right side corresponds to the case of Eq. (|T3|) . and "-" 
for Eq.flUD. 



Using results of papers [ill, ll2|, |l3j, Q, Q, the frequency of the breathing mode can be 
written as 

u 2 b = P/27rpRl (15) 
the frequency of the quadrupole surface modes mode as 

u 2 s = \2a/pRl (16) 

and the effective mass of the bubble as 

M = A7rpR 3 . (17) 

Here f3 is a constant dependant on the properties of the liquid, a is the surface tension, and 
p is liquid density. To estimate values of neutron cross sections on bubbles in liquid helium 
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FIG. 3: Inelastic total up-scattering cross sections of neutrons with energy E n = 100; neV as a 
function of a bubble radius in liquid helium for breathing (dashed line) and surface (solid line) 
modes, (n = 1) 



we use [Uj j3 = 1.02 -10" 1 eV nm~ 2 , p = 0.145 g cm~ 3 and a = 2.25 • 10 -3 eV nm~ 2 . In that 
case a bubble with a radius Rq = 1.75 nm has frequencies of oscillations fiu^ = 3.8 • 10~ 5 eV 
and huj s = 5 ■ 10~ 5 eV. Since these frequencies have the same order of magnitude, the 
main difference between total cross sections for the breathing and surface modes comes 
from the $- function in Eqs. ( jTTT) and ( |T2l) . Therefore, neutron cross sections have different 
energy dependence for breathing and surface modes(see figure ([3]); they also have different 
dependency on the bubble radius (see figure (j4])). 

Since the loss (or gain) of energy by neutron at each interaction is equal to hut, and 

—3/2 

noj oc Rq , this transferred energy could be adjusted by changing of the radius of the 
bubble. Due to theoretically understandable nature of neutron-bubble interactions, very 
cold and ultra cold neutrons could be to be used to study bubble formations and dynamics 
in liquids (for example, bubbles in liquid helium). However, the possible and even more 
important application of this process is the opportunity of the neutron moderation on the 
bubble gases. Inelastic scattering of neutrons on bubbles could slow down thermal neutrons 
to very cold or ultra cold energies by a small number of interactions without significant loss of 
the neutron intensity. It could increase the the rate of the production of very cold and ultra 



cold neutrons by orders of magnitudes. Moreover, using electron bubbles in liquid 



lie 



see, for example papers la, ll7|] and references therein) one can manipulate them 18 



mm 

] by 



S 



a(b) 

3.5; 

3: 

2.5- 
2 - 
1.5- 

1 , 
0.5- i 



/ 



/ 



40000 50000 60000 70000 80000 9000010000 



g(neV) 



FIG. 4: Inelastic total down-scattering cross sections on the bubble of radius of 1.75; nm in liquid 
helium as a function of neutron energy for breathing (dashed line) and surface (solid line) modes. 
(n = l) 

applying the proper electric and magnetic fields, which would lead to the opportunity 
to modulate fluxes of very cold and ultra cold neutrons. At the current stage, however, it is 
impossible to say how feasible these applications could be. They required more experimental 
studies, which are under consideration [19|. 
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